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Abstract.  Because  of  their  importance  in  the  accurate  modeling  of  the  Earth's  orientation  in  space,  some  non-negligible  pre¬ 
dictable  effects  on  precession-nutation  are  investigated.  This  paper  considers  the  coupling  effects  between  the  axial  and  the 
equatorial  components  of  the  Earth’s  rotation  vector  in  the  dynamical  equations,  and  the  effects  of  the  second  order  lunisolar 
torque  due  to  the  Earth’s  zonal  deformations.  Firstly,  the  coupling  effects  are  shown  to  contribute  for  less  than  0.1  // as  and  are 
therefore  negligible.  Secondly,  we  demonstrate  that  the  0.7  mas  contribution  of  the  rotation  rate  variations  due  to  zonal  tides 
to  the  nutation  in  obliquity  deduced  by  Bretagnon  et  al.  (2000,  Proc.  IAU  Coll.,  180,  230;  2001,  Celest.  Mech.  Dyn.  Astr.,  80, 
177)  is  an  artefact  which  comes  from  an  incomplete  way  of  taking  into  account  the  effect  of  the  rotation  rate  variations.  The  net 
contribution  is  shown  to  be  negligibly  small.  Thirdly,  for  an  Earth  model  with  an  elastic  mantle  and  decoupled  liquid  core,  the 
contribution  of  the  second-order  lunisolar  torque  due  to  the  Earth’s  zonal  deformations  is  shown  to  be  207.9  // as  and  -9.7  //as 
on  the  18.6-year  nutation  respectively  in  longitude  and  in  obliquity,  and  a  correction  of  — 4925.9  //as/century  on  the  precession 
in  longitude. 
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1.  Introduction 

Very  long  baseline  radio  interferometry  (VLBI)  observations 
provide  on  a  regular  basis  “celestial  pole  offsets”  that  con¬ 
tain  errors  in  the  models  for  the  position  of  the  celes¬ 
tial  pole  in  the  celestial  reference  system.  The  most  accu¬ 
rate  available  precession-nutation  model  is  the  IAU  2000A 
model  recommended  by  the  International  Astronomical  Union 
(IAU  Resolution  B1.8)  adopted  at  the  IAU  24th  General 
Assembly  (2000)  and  implemented  by  the  IERS  beginning  on 
1  January  2003. 

Comparison  of  this  new  model  against  VLBI  observations 
shows  differences  of  the  order  of  200  //as.  These  differences 
are  the  consequence  of  various  influences  as  well  as  geophysi¬ 
cal  processes  that  are  summarized  below.  The  main  oscillation 
appearing  in  these  residuals  comes  from  the  Free  Core  Nutation 
(FCN).  It  gives  rise  to  a  430-day  periodic  term,  the  amplitude 
of  which  varies  strongly  and  is  therefore  not  predictable  but  can 
be  described  using  an  empirical  model  (Herring  et  al.  2002).  As 
the  nutation  modeling  contains  only  permanent  predictable  ef¬ 
fects,  the  FCN  is  not  included  in  IAU  2000A  and  therefore  ap¬ 
pears  in  the  residuals.  Other  unpredictable  effects  come  from 
the  variable  atmospheric  and  oceanic  forcing  which  induces 


time  varying  nutation  amplitudes  of  about  80  //as  for  the  an¬ 
nual  term  (Bizouard  et  al.  1998).  Some  theoretical  predictable 
effects  are  taken  into  account  in  the  IAU  2000A  precession- 
nutation  model  whereas  they  were  not  included  in  the  previous 
models.  These  effects  include: 

-  the  couplings  between  axial  and  equatorial  components  ap¬ 
pearing  when  developing  the  dynamical  equations  of  the 
Earth  rotation  at  the  second  order; 

-  the  effects  of  variations  of  the  dynamical  ellipticity,  giving 
rise  to  changes  in  the  lunisolar  torque  exerted  on  the  Earth 
and  therefore  on  precession-nutation. 

Further  discussions  of  these  effects  are  still  necessary  to  clarify 
how  they  can  be  related  to  previous  studies.  Coupling  effects 
between  the  Earth’s  rotation  rate  and  precession-nutation  were 
neglected  because  their  contributions  are  expected  to  be  negli¬ 
gible  compared  to  the  accuracy  of  the  observations.  However, 
recent  studies  concluded  that  variations  in  Earth  rotation  due 
to  zonal  deformations  were  responsible  for  contributions  of  the 
order  of  712  //as  to  the  18.6-year  nutation  in  obliquity  and  a 
change  of  the  precession  in  longitude  of  4  mas  per  century 
(Bretagnon  et  al.  2000,  2001). 
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The  deformations  due  to  zonal  tides  have  also  direct  ef¬ 
fects  on  precession  and  nutation.  The  effect  on  nutation  was 
first  studied  by  Souchay  &  Folgueira  (2000).  It  was  concluded 
that,  for  a  basic  Earth  model,  this  should  modify  the  amplitude 
of  the  18.6-year  nutation  in  longitude  by  168  /./as,  perturba¬ 
tions  on  other  nutation  terms  all  being  less  than  1 0  /./as.  Another 
work  by  Krasinsky  (1999)  showed  larger  corrections  based  on 
a  computation  of  the  gravitational  torque  exerted  on  the  Earth. 
The  largest  correction  provided  is  5.11  mas  on  the  13.66-day 
term,  whereas  the  semi-annual  term  and  the  18.6-year  term 
are  changed  respectively  by  2.38  mas  and  1.07  mas.  Note  that 
VLBI  observations  do  not  show  such  large  discrepancies  in  the 
nutation  model. 

Considering  such  a  large  disparity  in  the  amplitudes  of  the 
contributions  to  nutation  due  to  zonal  tides  that  are  reported  in 
various  studies,  a  rigorous  computation  of  these  effects  is  nec¬ 
essary.  In  this  paper,  we  compute  the  effects  due  to  the  torque 
produced  by  the  action  of  the  tesseral  potential  on  the  defor¬ 
mations  produced  by  the  zonal  potential  together  with  the  ef¬ 
fects  on  the  rotation  rate  produced  by  these  deformations.  In 
Sect.  2  we  recall  the  expressions  of  the  Euler-Liouville  equa¬ 
tions  developed  up  to  the  second  order,  emphasing  the  cou¬ 
plings  between  variations  of  the  axial  component  and  the  equa¬ 
torial  component  of  the  instantaneous  rotation  vector.  In  Sect.  3 
we  introduce  the  second  order  part  of  the  lunisolar  torque  com¬ 
ing  from  Earth’s  zonal  deformations,  and  in  Sect.  4,  we  com¬ 
pute  all  the  predictable  effects  due  to  zonal  tides  on  the  Earth’s 
precession-nutation. 


2.  Dynamical  equations  of  the  Earth’s  rotation 
up  to  the  second  order 

We  consider  a  deformable  Earth  with  fluid  core.  In  a  rotating 
frame  with  its  axes  oriented  towards  the  Earth’s  mean  axes  of 
inertia,  the  angular  momentum  equation  is: 


d  H 

dr 


+  Cl)  X  H  —  r 


(1) 


where  H  is  the  whole  Earth’s  angular  momentum  vector,  cor¬ 
responding  to  the  product  of  the  inertia  tensor  7: 
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wherein  second-order  quantities  c,;  characterize  the  effects  of 
the  deformations,  with  the  Earth’s  instantaneous  vector  of  rota¬ 
tion  a>: 


a>  —  Q 


m  i 
m2 

1  +  m3 


(3) 


Q  being  the  mean  Earth’s  rotation  rate, 

r  is  the  tidal  torque  expressed  in  the  terrestrial  frame.  The 
equatorial  component  T  =  Ti  +  UN  of  the  torque  is  computed 
from  the  expression  P  in  the  celestial  frame  by  the  relation: 

T  =  TV1®  (4) 


and  the  Euler’s  kinematical  relations: 

8  +  TT  sin  8  =  -Lime1®  (5) 

<t>  +  T1  cos  8  =  0(1+ m3) 


relate  the  terrestrial  motion  of  the  Earth’s  axis  of  rotation 
m  —  m\  +  i/m  and  the  relative  variations  of  the  Earth’s  rotation 
rate  m3  to  the  Euler’s  angles  8,  'F  and  <1>  between  the  terrestrial 
frame  and  the  mean  ecliptic  frame  of  the  epoch  of  reference. 

Using  the  equations  of  Sasao  et  al.  (1980),  one  obtains  the 
following  expressions: 
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in  which  c  —  C13  +  ic23,  c'f  =  Cf,i3  +  icf,23,  0+  =  0(C  -  A) /A  and 
/Tf  =  0(Cf  -  Af)/Af,  quantities  subscripted  by  f  being  relevant 
to  the  core. 

The  previous  equations  reflect  two  different  couplings  be¬ 
tween  the  axial  and  the  equatorial  component  of  the  rotation 
vector  of  the  Earth.  One  coupling  comes  from  the  products 
m  x  m3  or  nif  x  m3  in  the  left  hand  side  of  (6).  Another  cou¬ 
pling  appears  through  the  sidereal  rotation  angle  depending 
upon  m 3,  in  Eqs.  (4)  and  (5)  in  the  opposite  sense  to  transform 
coordinates  between  the  terrestrial  to  the  celestial  frame. 

Variations  of  Euler’s  angles  are  related  to  the  instantaneous 
vector  of  rotation  by  Euler’s  kinematical  relations  (5).  One  gets 
the  nutation  angles  referred  to  the  ecliptic  of  epoch: 


A1//1  sin  ei  =  AT'  sin  8 
Aei  =  -A  8. 


(7) 


Note  that  the  periodic  part  of  the  nutation  angles  referred  to 
the  ecliptic  of  date  (At//  and  Ae)  is  the  same  as  the  one  of  Ai//\ 
and  Aei,  and  that  the  precession  is  expressed  along  the  fixed 
ecliptic.  The  contribution  to  the  precession  1//3  is  equivalent  to 
the  secular  term  of  Ai//\ . 


3.  Torque  including  Earth’s  zonal  deformations 

The  tidal  gravitational  torque  exerted  on  the  Earth  is  derived 
from  the  tesseral  part  of  the  lunisolar  potential  </>  in  the  celestial 
reference  frame  (Sasao  et  al.  1980): 

r  =  r"1^  +  r(2)  (8) 

=  -i Qr(C  -  A)cf>  -  i£22(c33  -  cu)<P 


where  the  superscripts  (1)  and  (2)  stand  for  the  first  and  second 
order  parts  of  the  torque  respectively.  The  second  order  part 
includes  the  deformations  C33  and  cn  of  the  Earth  along  the 
directions  of  the  mean  axes  of  inertia  (cf.  relation  (2)). 

The  complex- valued  quantity  </>  depends  on  the  right  ascen¬ 
sion  a  and  declination  5  of  the  perturbing  body: 


0  = 


3K2/j 
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sin  c5  cos  (/(cos  a  +  i  sin  a) 


(9) 
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Table  1.  Excess  of  LOD  (yrs)  due  to  zonal  tides  from  IERS  Conventions  2003,  relative  variations  of  the  Earth’s  rotation  rate  mi  (10  12  rad)  and 
C33/C  (p as)  from  Eqs.  (11)  and  (12). 
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in  which  K 2  is  the  geocentric  gravitational  constant,  p  is  the 
ratio  between  the  mass  of  the  disturbing  body  and  the  mass  of 
the  Earth  and  r  is  its  geocentric  distance. 

Assuming  as  generally  (see  e.g.  Melchior  1978;  Yoder  et  al. 
1981)  that  the  trace  of  the  inertia  tensor  is  constant,  one  has  in 
the  case  of  zonal  deformations: 

6{C  -  A)  =  C33  -  cii  =  ^c33.  (10) 

The  variations  of  the  dynamical  ellipticity  are  related  to  corre¬ 
sponding  changes  in  the  excess  of  length-of-day  (LOD)  and  to 
the  axial  component  m3  of  the  instantaneous  vector  of  rotation 
of  the  whole  Earth: 


LOD 


with  a  microarc  second  accuracy  requires  a  more  sophisticated 
Earth  model.  The  effect  of  the  Earth  model  appears  through  the 
zonal  response  coefficient  that  can  be  estimated  from  observed 
Earth  Orientation  Parameters  time  series  (Hefty  &  Capitaine 
1990).  The  model  of  zonal  variations  in  Earth’s  rotation  of 
Yoder  et  al.  (1981)  considered  an  equilibrium  ocean  tide  and 
was  adopted  in  the  IERS  Conventions  1996  (McCarthy  1996). 
A  more  refined  model  provided  by  Defraigne  &  Smits  (1999) 
was  adopted  in  the  IERS  Conventions  2003.  It  includes  a  de¬ 
coupled  core  and  an  inelastic  mantle,  and  ocean  effects  in¬ 
cluded  using  a  transfer  function  worked  out  by  Mathews  et  al. 
(2002). 

Our  computation  of  C33  uses  this  model  for  LOD  variations 
and  the  mean  moments  of  inertia  of  the  whole  Earth  provided 
by  the  IERS  Conventions  2003: 


where  the  superscript  z  indicates  that  these  rotation  rate  varia¬ 
tions  come  from  zonal  tides  and  LOD  =  86400  s.  According 
to  Mathews  et  al.  (2002)  and  Mathews  (2004),  the  increment  of 
inertia  of  the  whole  Earth  produced  by  the  deformations  due  to 
zonal  tides  is: 


C33  =  -CeSm'3  (12) 

r  — 

eff  “  1  -  yCf/KC 

y  and  k  being  undimensional  parameters  representing  the  de- 
formabilities  of  the  mantle  and  the  core-mantle  boundary  un¬ 
der  degree  2  tidal  forcing  (see  e.g.  Sasao  et  al.  1980;  Mathews 
et  al.  1991)  and  Cf  the  axial  moment  of  inertia  of  the  core. 

The  second  order  torque  due  to  the  zonal  deformations  is 
finally  expressed  using  quantity  C33  or  equivalently  m\: 


„,(2)  .9  K2n  .  .  ..  .  .  . 

I  K  ’  =  -1  -  -  C33  sin  0  cos  o(cos  a  +  1  sin  a) 


2r3 
.9 gVeff 
2  r3 


(13) 


sin  6  cos  d(cos  a  +  i  sin  a). 


Direct  computation  of  the  changes  in  the  dynamical  ellipticity 
can  be  based  on  expression  for  C33  corresponding  to  a  given 
Earth  model.  For  example  Souchay  &  Folgueira  (2000)  used 
a  model  provided  by  Melchior  (1978)  on  the  basis  of  a  sim¬ 
plified  elastic  Earth  model.  Computing  the  effect  on  nutation 


C  =  8.037  x  1037  kg  nr  (14) 

A  =  8.011  x  1037  kg  m2 
Cm  =  7.040  x  1037  kg  m2. 

Using  the  compliances  from  Mathews  et  al.  (2002)  estimates, 
one  gets: 

Ceff  =  9.213  xlO37  kg  nr.  (15) 

Table  1  gives  the  excess  of  the  length-of-day  due  to  tidal  defor¬ 
mations  provided  by  the  IERS  Conventions  2003  together  with 
the  variations  mi.  and  the  relative  variations  of  C  from  Eqs.  (11) 
and  (12). 


4.  Solution  of  the  dynamical  equations 

4.1.  Couplings  with  the  rotation  rate  m3 

All  the  terms  showing  products  m  x  m3  or  mt  x  m3  in  the  left- 
hand  side  of  Eq.  (6)  express  coupling  between  the  motion  of 
the  axis  of  rotation  and  the  Earth’s  rotation  rate  (note  that  the 
use  of  a  decoupled  core  involves  the  relation  nify  -  -m3).  The 
largest  contribution  of  these  terms  to  nutation  has  been  evalu¬ 
ated  as  being  lower  than  0.1  pas  for  the  18.6-year  term.  This 
confirms  that  this  coupling  effect  is  negligible  and  it  will  not  be 
considered  further  in  the  following. 
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4.2.  The  coupling  with  the  sidereal  rotation  angle  <P 


The  effect  of  taking  into  account  the  Earth’s  rotation  irregular¬ 
ities  in  the  sidereal  rotation  angle  <5  in  computing  the  solution 
for  precession-nutation  has  to  be  investigated. 

Bretagnon  et  al.  (2000,  2001)  constructed  non-rigid  nu¬ 
tation  series  by  applying  the  transfer  function  MHB  2000 
(Mathews  et  al.  2002)  to  the  SMART97  rigid  nutation  series 
(Bretagnon  et  al.  1998).  Given  that  the  transfer  function  of 
MHB  2000  does  not  take  into  account  the  variations  in  the 
Earth’s  rotation  rate,  these  effects  were  added  to  the  third  com¬ 
ponent  of  the  Earth’s  rotation  vector  and  the  Euler’s  angles 
were  recomputed  from  Euler’s  kinematical  relationships.  This 
method  gave  rise  in  particular  to  a  18.6-year  term  in  obliquity 
with  an  amplitude  reaching  712  //as. 

We  investigate  this  effect  considering  the  equations  for  a 
rigid  ellipsoidal  Earth  (the  quantity  c  is  null)  together  with  vari¬ 
ations  in  Earth’s  rotation  rate.  Equation  (6)  is  reduced  to: 


m  —  i  crrm 


r 

AO 


and  therefore: 


(16) 


m(t)  =  Kte1<Tr‘  +  ei<Tri'  f  —  e^V^df  (17) 

J  AO 

where  Kr  is  a  constant  depending  on  initial  values  of  the  prob¬ 
lem. 

The  first  term  of  solution  (17)  is  the  free  motion  of  the  rota¬ 
tion  axis  of  the  Earth,  known  as  the  Euler  motion.  The  second 
term  of  Eq.  (17)  corresponds  to  the  forced  motion  of  the  rota¬ 
tion  axis.  Variations  of  the  Euler’s  angles  are  therefore: 


e  +  i'f'sine  =  -fWV’*  f  —  e^V^'df 

J  AO 


<f>  +  vFcos0  =  0(1+ m3). 


(18) 


Using  an  iterative  process,  the  integration  of  the  system  (18) 
converges  rapidly  after  a  few  iterations.  Note  that  the  side¬ 
real  rotation  angle  <1>  is  affected  by  the  irregularities  >m  in  the 
Earth’s  rotation  rate,  according  to  the  second  equation  of  the 
system.  In  the  first  equation  of  Eq.  (18),  the  quantity  <5  appears 
twice  inside  and  outside  the  time  integral  in  opposite  sense. 
First,  it  is  used  to  transform  the  torque  T'  from  the  celestial 
frame  to  the  terrestrial  frame  (rotation  of  -<f>).  Second,  it  is 
used  in  Euler’s  kinematical  relations  to  transform  the  rotation 
vector  expressed  in  the  terrestrial  frame  into  the  Euler’s  an¬ 
gles  with  respect  to  the  space-fixed  reference  frame  (rotation 
of+O). 

By  a  numerical  evaluation,  one  can  see  that  omitting  the 
variations  in  the  Earth’s  rotation  rate  in  the  two  rotations  of  an¬ 
gle  plus  or  minus  <1>,  that  is  to  say  writing  <i>  +  V  cos  6  =  Q. 
instead  of  <j>  +  vt'cos  6  =  Q(1  +  m3),  leads  to  an  error  below 
0. 1  //as  on  6  and  ML  The  effect  of  the  former  rotation  is  actu¬ 
ally  mostly  cancelled  out  by  the  latter  so  that  the  net  effect  is 
negligibly  small. 

However,  omitting  to  take  into  account  the  Earth’s  rotation 
rate  irregularities  in  one  of  the  rotations  of  plus  or  minus  <1> 
leads  to  large  terms  which  could  be  interpreted  as  being  cou¬ 
pling  effects.  For  example,  if  the  computation  of  the  torque  in 


Table  2.  Effect  of  zonal  variations  in  Earth's  rotation:  artefacts  in¬ 
duced  by  omitting  the  variations  when  transforming  the  torque  from 
the  celestial  frame  to  the  terrestrial  frame  compared  to  the  values  of 
Bretagnon  et  al.  (2001  j,  and  net  contribution. 


Artefact 

Net  effect 

Ai// 

Ae 

A/A 

Ae 

//as 

//as 

//as 

//as 

1  r  f  d  n 

sin 

cos 

sin 

cos 

0000  1 

This  study 

4 

-729 

1 

-1 

Bretagnon  et  al.  (2001) 

4 

-712 

0  0  0  0  2 

This  study 

71 

19 

0 

0 

Bretagnon  et  al.  (2001 ) 

69 

19 

0  0  2  -2  3 

This  study 

9 

-3 

0 

0 

Bretagnon  et  al.  (2001 ) 

- 

- 

002-2  1 

This  study 

-8 

3 

0 

0 

Bretagnon  et  al.  (2001 ) 

- 

- 

Precession  in  longitude 

//  as/c 

This  study 

4497 

Bretagnon  et  al.  (2001) 

4224 

the  terrestrial  frame  neglects  these  irregularities,  one  gets  a  sub¬ 
stantial  contribution  in  the  in-phase  1 8 . 6-year  term  of  the  obliq¬ 
uity  (-729  //as)  and  on  the  9.3-year  in-phase  term  of  the  longi¬ 
tude  (7 1  //as).  This  1 8.6-year  oscillation  of  the  angle  9  is  due  to 
a  coupling  between  the  precession  in  longitude  T  =  To  x  t  and 
the  18.6-year  term  in  the  variations  of  the  Earth’s  rotation  rate 
of  amplitude  a\s.6- 

9i&.6  =  'Po  sin  8  6  cos<X|8  6/.  (19) 

°"l8.6 

The  precession  in  longitude  is  also  affected  by  a  spurious  ef¬ 
fect  of  4  mas/c.  This  explains  the  erroneous  values  provided  in 
the  above-mentioned  studies  that  followed  such  an  approach. 
These  “artefacts"  are  displayed  in  Table  2  together  with  the 
terms  provided  by  Bretagnon  et  al.  (2000,  2001)  and  the  net 
effect. 

4.3.  Effect  of  the  second  order  zonal  tides-induced 
torque 

The  computation  of  the  effects  of  the  second  order  torque  in¬ 
duced  by  zonal  tides  on  nutation  is  done  using  Eq.  (6)  in  the 
frequency  domain  once  the  torque  is  computed  according  to 
Eq.  (13).  The  calculation  procedure  of  the  lunisolar  torque 
made  in  the  present  paper  is  based  on  the  lunar  theory  ELP2000 
(Chapront-Touze  &  Chapront  1983)  and  the  solar  system  semi- 
analytical  solution  VSOP87  (Bretagnon  &  Francou  1988). 
Poisson  series  manipulations  have  been  processed  using  the 
GREGOIRE  software  package  developed  by  Chapront  (2003). 
We  use  standard  values  (1ERS  Conventions  2003)  for  other  or¬ 
bital  parameters  involved  in  the  computation  of  the  torque. 


S.  Lambert  and  N.  Capitaine:  Zonal  deformations  and  precession-nutation 


259 


The  computed  variations  of  the  nutations  angles  are  the  fol¬ 
lowing: 

(i)  Secular  term.  The  constant  term  of  the  torque  leads  to  a  sec¬ 
ular  variation  in  longitude  of  -4926  //as/c.  Note  that  such 
a  contribution  would  be  automatically  included  in  the  ob¬ 
served  value  for  the  precession.  The  obliquity  is  not  af¬ 
fected. 

(ii)  Periodical  terms.  The  main  second-order  contribution  to 
the  periodic  terms  comes  from  the  18.6-year  term  of  the 
torque  which  is  a  coupling  between  the  18.6-year  vari¬ 
ation  of  the  dynamical  ellipticity  and  the  constant  term 
of  the  lunisolar  potential.  Other  contributions  are  below 
5  //as.  It  is  interesting  to  note  that,  according  to  Table  1,  the 
1 8.6-year  term  is  in  the  opposite  sense  compared  to  the  oth¬ 
ers.  Namely,  it  is  negative  and  decreases  the  excess  of  LOD 
(or  equivalently,  decreases  the  difference  between  axial  and 
equatorial  moment  of  inertia  so  that  the  Earth  becomes  less 
elliptic).  This  means  that  the  effect  of  the  18.6-year  vari¬ 
ation  is  to  concentrate  the  masses  of  the  Earth  around  the 
rotation  axis,  so  that  the  Earth  rotation  speed  increases  to 
satisfy  the  angular  momentum  conservation  law,  contrar- 
ily  to  the  other  contributions  which  are  positive  and  make 
the  Earth  oblateness  increase.  A  consequence  of  the  oblate¬ 
ness  decrease  is  that  the  torque  exerted  on  the  equatorial 
bulge  is  lowered  and  the  18.6-year  nutational  response  to 
this  excitation  will  consequently  be  lowered.  The  18.6-year 
nutation  coefficient  should  be  decreased  in  absolute  value. 
In  fact,  this  contribution  reaches  208  //as  in  the  nutation  in 
longitude  whereas  it  produces  only  a  small  effect  of  almost 
-10  //as  on  the  nutation  in  obliquity.  Since  the  18. 6-year 
nutation  amplitude  in  longitude  is  negative  (-17"),  such  a 
result  is  in  the  right  direction. 

Table  3  summarizes  the  results  and  provides  a  comparison  be¬ 
tween  our  study  and  other  works.  The  difference  of  the  order 
of  40  //as  for  the  18.6-year  nutation  in  longitude  with  respect  to 
Souchay  &  Folgueira  (2000)  is  due  to  the  more  realistic  Earth 
model  considered  in  this  study.  The  differences  with  respect  to 
Mathews  et  al.  (2002)  need  some  additional  comments. 

Mathews  et  al.  (2002)  investigated  a  refined  Earth  model 
(elastic  mantle  and  decoupled  liquid  core)  to  evaluate  the  in¬ 
fluence  of  the  so-called  "non-linear”  terms  in  the  dynamical 
equations  in  order  to  remove  them  from  the  observations  before 
fitting  the  geophysical  parameters  used  in  the  transfer  function. 
This  gave  a  contribution  to  the  18.6-year  term  of  94  //as  and 
-29  //as  respectively  in  longitude  and  in  obliquity,  including 
the  second  order  torque  due  to  both  zonal  and  sectorial  tides. 
Note  that  the  differences  between  the  values  obtained  in  this 
study  and  those  of  Mathews  et  al.  (2002)  appearing  in  Table  3 
are  in  fact  considerably  reduced  with  respect  to  revised  values 
obtained  by  Mathews  (2004).  Excluding  the  sectorial  part,  the 
revised  contributions  in  the  18.6-year  nutation  are  194  //as  and 
-10  //as  respectively  in  longitude  and  in  obliquity  (Mathews 
2004).  The  remaining  difference  with  respect  to  our  values  is 
due  to  anelasticity  and  ocean  tide  effects  taken  into  account 
in  the  MHB  work.  It  should  moreover  be  noted  that  since  this 
work  has  been  carried  out,  Mathews  (2003)  pointed  out  that  the 


Table  3.  Effects  of  the  second  order  contribution  of  the  tidal  poten¬ 
tial  (zonal  tides  contribution)  on  nutation  angles  (//as)  and  precession 
(// as/c)  from  different  studies.  MHB  2000Z:  values  of  Mathews  et  al. 
(2002)  including  only  the  zonal  tides  contribution;  SF99:  Souchay  & 
Folgueira  (2000). 


Ai/'  Ae 

//  as  //as 


/ 

r 

FDD. 

sin 

COS 

sin 

cos 

0 

0 

0  0  1 

MHB  2000Z 

*101.3 

- 

- 

*-8.9 

SF99 

168.0 

- 

- 

-9.0 

This  study 

207.9 

2.4 

0.2 

-9.7 

0 

0 

0  0  2 

MHB  2000Z 

- 

- 

- 

- 

SF99 

4.5 

- 

- 

-2.8 

This  study 

5.0 

0.1 

0.0 

-3.2 

0 

0 

2-2  2 

MHB  2000Z 

- 

- 

- 

- 

SF99 

3.2 

- 

- 

0.6 

This  study 

4.0 

0.0 

0.0 

0.4 

Precession  in  longitude 

//as/c 

MHB  2000Z 

- 

SF99 

- 

This  study 

^1925.9 

*Note  that  revised  values  by  Mathews  (2004)  are  194  //as 
and  -10  //as  for  the  18.6-year  nutation  in  longitude  and  obliquity,  re¬ 
spectively. 


above  computed  effects  could  be  cancelled  out  by  the  sectorial 
and  zonal  parts  of  the  potential  acting  on  the  deformations  due 
to  the  tesseral  potential.  This  point  needs  further  investigation. 


5.  Discussion  and  conclusion 

Previous  studies  (Bretagnon  et  al.  2000,  2001)  pointed  out  that 
the  Earth’s  rotation  rate  variations  due  to  zonal  tides  could  have 
a  noticeable  effect  on  precession-nutation,  with  an  amplitude 
of  the  order  of  700  //as  in  the  18.6-year  nutation  in  obliquity, 
although  such  an  effect  was  not  detected  in  the  observations. 

The  present  paper  has  shown  that  this  effect  is  in  fact  can¬ 
celled  out  by  taking  into  account  the  variations  of  the  rotation 
rate  in  the  expression  of  the  torque  as  seen  from  the  terrestrial 
frame  used  in  the  dynamical  equations.  This  point  resolves  a 
serious  dilemma  of  several  years  concerning  the  magnitude  of 
these  effects.  Moreover,  this  term  does  not  contain  the  effects 
induced  by  the  variations  of  the  dynamical  ellipticity,  although 
it  may  sometimes  have  been  understood  as  being  so. 

We  also  recomputed  the  effects  due  to  the  coupling  be¬ 
tween  axial  and  equatorial  components  of  the  rotation  vector 
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in  the  dynamical  equations  of  the  Earth  rotation  and  we  took 
into  account  the  second-order  contribution  of  the  torque  due  to 
zonal  deformations  of  the  Earth’s  shape.  These  results  are  ob¬ 
tained  by  solving  the  dynamical  Eqs.  (6)  using  the  model  for 
variations  in  the  rotation  rate  due  to  zonal  tides  of  the  IERS 
Conventions  2003.  For  an  elastic  Earth  with  a  decoupled  liquid 
core,  we  have  concluded  that: 

•  the  coupling  effects  between  axial  and  equatorial  compo¬ 
nents  of  the  rotation  vector  do  not  bring  any  contribution 
larger  than  0. 1  /./as  and  are  therefore  negligible; 

•  the  second-order  torque  induced  by  Earth’s  zonal  deforma¬ 
tions  produces  larger  effects.  The  effect  on  the  18.6-year 
nutation  is  of  207.9  //as  in  longitude  and  -9.7  //as  in 
obliquity.  Other  contributions  are  about  a  few  microarc- 
seconds.  The  correction  to  the  precession  in  longitude  is 
-4925.9  //as/c. 

Note  that  for  comparison  with  VLBI  observations,  the  other 
second  order  contributions  of  the  torque  on  precession-nutation 
as  mentioned  in  Mathews  (2003)  should  be  considered. 
Although  the  effect  of  zonal  tides  on  nutations  is  already  in¬ 
cluded  in  MHB  2000  and  should  be  completed  by  other  second 
order  effects  (not  yet  published),  this  computation  gives  a  more 
realistic  model  than  that  of  Souchay  &  Folgueira  (2000)  as  well 
as  a  result  independent  of  MHB  2000  and  confirms  an  inaccu¬ 
racy  in  the  MHB  2000  tables. 
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